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Abstract This paper deals with the kinetics of the color 
changes of cellulose during heat treatment. The color of 
cellulose heated at 90-180°C was measured by a spectro¬ 
photometer and expressed by CIELAB color parameters. 
The values of L* decreased and those of a*, b* and AE* ab 
increased at all the treatment temperatures. Several kinetic 
models, namely, the zero-order, first-order, second-order 
and autocatalytic model, were applied to the changes in the 
color values. Furthermore, the results of kinetic analysis 
using the best-fit model were compared to the results 
obtained from conventional kinetic models. It was sug¬ 
gested that the kinetic analysis using the best-fit model was 
the better way to accurately predict color changes during 
heat treatment. The values of apparent activation energy 
calculated from the changes of L*, a*, b* and AE* ab were 
125, 124, 118 and 120 kJ/mol, respectively. These values 
were similar to the reported values calculated from other 
chemical or mechanical properties. 
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Introduction 

Cellulose is the major component of ligno-cellulosic 
materials such as paper, textile and wood, and the 
chemical and mechanical degradation of cellulose has a 
great influence on the degradation of ligno-cellulosic 
materials. Kinetic study has been of great importance not 
only for understanding the degradation mechanism but 
also for predicting the degradation of chemical and 
mechanical properties. Numerous kinetic studies on 
cellulose have been reported using thermally and hygro- 
thermally treated cellulose, mainly focusing on the 
degradation of electrical insulation paper in electrical 
power transformers [1, 2] and the degradation of acidic 
paper, which has been widely used from the 18th century 
and has caused catastrophic degradation of books and 
documents [3-5]. 

The degradation of cellulose has been widely elucidated 
by weight decrease, degree of polymerization (DP), 
mechanical properties, crystallinity and morphological 
characteristics. In addition, the color difference and the 
yellowing index are also used [6, 7]. In the field of paper 
and textile science, the change in color of cellulose is 
important not only from an aesthetic point of view but also 
as a possible index of the degradation of materials [8]. 
Several investigations have reported a correlation between 
color parameters and either DP or mechanical properties 
[9-11]. Color measurement is non-destructive and only 
requires a small area, unlike the measurements of chemical 
and mechanical properties. Therefore, color measurement 
is also suitable when the aesthetic aspects of measuring 
objects are important and destructive tests are unaccept¬ 
able, for example, for the investigation of culturally sig¬ 
nificant properties or the quality monitoring in industrial 
processing. 
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Table 1 The basic information of the cellulose sample and the average initial color values with their standard deviations in parentheses 


Trade name 

Ash content (%) 

Basis weight (g/m") 

Degree of polymerization 

T 0 * «o* ^o* 

Whatman No.l 

<0.06 a 

87 a 

500-5000 b 

98.11 (0.01) -0.02 (0.01) 2.31 (0.06) 


a Based on the company’s information 
b Based on the literature [15-18] 


The color of cellulose turns brown or yellow during heat 
treatment [11, 12]. However, the mechanism of this change 
has not been clarified yet and only a few reports were found 
in kinetic analysis based on the color changes of cellulose 

[13]. 

This study examines the color changes of cellulose 
during heat treatment to elucidate and predict desirable/ 
undesirable color changes of ligno-cellulosic materials. For 
an accurate analysis, we have applied several kinetic 
models and a best-fit model employing a kinetic approach 
based on the time-temperature superposition (TTSP) 
method. 

Materials and methods 

Materials 

Whatman No.l cellulose filter paper was used as a cellulose 
sample. The basic information of the samples is summa¬ 
rized in Table 1. Specimens of approximately 20 x 40 mm 
were dried in an air-circulating oven and then in a vacuum 
oven at 60°C for 24 h each. The low drying temperature 
was set to avoid undesirable reactions of cellulose. Each 
dried specimen was heated in an air-circulating oven at 4 
temperature levels from 90 to 180°C for a duration ranging 
from an hour to approximately 1.8 years. Table 2 shows 
the treatment times of the specimens at each treatment 
temperature. They were planned under the assumption that 
a 10°C increase of the reaction temperature halves the 
reaction time to obtain a similar intensity of reaction among 
the treatment of different temperatures [14]. The treatment 
at 90°C is now in process with the planned treatment time 
shown in parentheses in Table 2. 

Color measurement 

The color of the specimens was measured with a spectro¬ 
photometer (KONICA MINOLTA CM-2600d) under a 
D65 light source and an observed angle of 10°. The sensor 
head of the spectrophotometer was 8 mm in diameter. The 
size of paper specimens was sufficiently large relative to 
the sensor head. The CIELAB color parameters (L*, a* and 
b*) were used to express the color changes. The differences 
in parameters AL*, A a* and Ab* and the total color 


Table 2 Treatment temperatures and times of cellulose samples 


90°C 

120°C 

150°C 

180°C 

Treatment time (h) 

64 

8 

1 

768 






16 

2 

1560 

192 

24 

3 

2560 

320 

40 


3584 

448 

56 

7 

6144 


96 


9216 




12288 

1536 

192 

24 

15580 

3096 

384 

48 

(24576) 

5450 

672 


(43008) 

7680 

960 

120 

(61440) 


1280 





240 



2304 

288 



2768 

336 

The treatment at 90°C is now in process with the planned treatment 


time in parentheses 

differences AE* ab were calculated using the following 
formulas: [19, 20] 

AV = L* - L* 

Aa*=a*-al (1) 

Ab* =b* -b* 0 

AE* b = \JAL* 2 -j- Aa * 2 -j- Ab * 2 (2) 

where L* is the lightness and a* and b* are the color 
coordinates under any testing condition, respectively, and 
L 0 *, a Q * an d b$* are the corresponding reference values 
obtained as the average of untreated specimens shown in 
Table 1, respectively. The color values of 3 specimens 
were averaged and the standard deviations were also cal¬ 
culated. The average of the coefficients of variation for 
each treatment condition were 0.23, 2.70, 1.38 and 1.89% 
for L*, a*, b* and AE* ab , respectively. Since the variation 
in the color values was negligibly small, the error bars 
displaying standard deviations were not drawn in the fol¬ 
lowing figures. 

The specimen preparation and the color measurement 
were always performed using rubber gloves. 
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Weight reduction 

The sample amount for each condition was approximately 
3 g, which was measured before and after heat treatment 
under the dry condition. The weight reductions were cal¬ 
culated based on the following equation: 

Wn - W t 

Weight reduction = 100 x- '-(%) (3) 

W 0 

where W 0 and W t are the dry weight of cellulose before 
treatment and at treatment time t , respectively. 

Kinetic analysis for color changes of cellulose 

Numerous reports have discussed the applicability of 
kinetic models, especially with respect to the degree of 
polymerization and some mechanical properties [21, 22]. 
However, there have been few reports dealing with the 
kinetic analysis of color properties [13]. 

The rate of change of a material property P at time t can 
be represented by: 

d P 

— = ±kP n (4) 

d t 

where n is the order of reaction. 

By integrating Eq. (4), the zero-order (Eq. 5), first-order 
(Eq. 6) and second-order (Eq. 7) kinetic models can be 
derived as: 


Pq dz kt 

(5) 

P 0 exp (±kt) 

(6) 

Po 

1 dz ktPo 

(7) 


where P is the color value at a treatment time t , P 0 the 
initial value of color, and k the rate constant of the 
chemical reaction. Many cellulose degradation kinetics 
have been based on these reaction models. Zervos et al. 
[13] reported that a kinetic model based on autocatalysis 
(Eq. 8) fitted well the changes of L* and Z?* during the heat 
treatment in sealed vessels: 

P = P 0 ± c( 2 kr - 1) (8) 

where C is constant. 

In general, it is known that a chemical reaction can be 
described by the Arrhenius equation: 

k = A exp ( 9 ) 

where A is the frequency factor, E a the apparent activation 
energy, R the gas constant, and T the absolute temperature 
of the reaction. The apparent activation energy is then 
obtained from the slope of the Arrhenius plot, which is the 
logarithm of k versus the reciprocal of T [23]. 


To determine the activation energy more accurately, a 
TTSP method has been proposed [24-27]. An Arrhenius 
approach applying the TTSP method is briefly described 
below. When the measured material property is plotted as a 
function of logarithmic treatment time at different tem¬ 
peratures, the shift distance used to superpose the plotted 
curves to the curves at a fixed reference temperature T ref 
along the logarithmic time axis is called the time- 
temperature shift factor a T , and is given as: 

cij = — (10) 

^ref 

where t ref is the test time at T ref , and t T is the time required 
to give the same response at the test temperature 
T. Combining Eqs. (9) and (10) gives: 


<2t = exp 






where both T and T ref are absolute temperatures. Plotting 
In (a T ) versus 1/7" is another way to calculate E a . 

Data analysis was performed using the FlexPro 8 soft¬ 
ware package (Weisang GmbH, Germany). All parameters 
of the model functions and the values of a T were estimated 
using a non-linear iterative curve-fitting method. The 
coefficient of determination (R ) and the root mean square 
error (RMSE) were used as the criteria to express the fitting 
of the tested models to the experimental data: 


*2 EL(*w ) 2 

R = —* -" 


EL (y { - ry 


RMSE = 


EL (Yi - Yi) 

N 


( 12 ) 

(13) 


A 

where Y t and Y t are the measured and modeled values, 
respectively, Y is the average of the measured data and N is 
the total number of observations. The higher the value of 
R ~ and the lower the value of RMSE, the better the model 
fits the data. 


Results and discussion 

Color changes of cellulose during heat treatment 

As shown in the representative specimen presented in 
Fig. 1, the color of cellulose turned brown and dark with 
increasing treatment time. Figure 2 shows the changes of 
the color parameters according to the duration and tem¬ 
perature of treatment. The decrease in L* at all tempera¬ 
tures indicates that the specimens became darker with 
increasing treatment time. The increases in a* and b* 
indicate that the specimens became more red and yellow. 
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Fig. 1 Cellulose specimens 
treated at 180°C with the 
treatment time (h) 



Fig. 2 Changes in L*, < 2 *, b*, 
and AE* ab during heat treatment 
as a function of time 




log (Treatment time) [h] log (Treatment time) [h] 



log (Treatment time) [h] 



log (Treatment time) [h] 


These results corresponded to the increase in AE* ab . 
Similar behavior irrespective of treatment temperatures 
suggested that the same reaction processes might cause the 
color changes, irrespective of temperatures. 

Kinetic analysis of the color changes of cellulose 

Several kinetic models were applied to the experimental 
data. In addition, in an analysis independent of these 
models, the function giving the best fit was selected. Since 
a similar behavior was observed irrespective of treatment 
temperatures, the data from the 180°C treatment were used 
to determine the best-fit model. 

Judging from the values of R and RMSE, the auto- 
catalytic model fit the color changes well in all color 
parameters, while the zero-, first- and second-order kinetic 
models did not fit the color parameters well except in the 
case of L*. 

As the function giving the best fit for all color param¬ 
eters, the NIST Hahnl model was selected: 


p Cp + Ctr + Cz^ + Cjf 3 

1 + c 4 t + C 5 t 2 + c 6 fi [ ’ 

where Q (i = 0, ..., 6) is constant. However, it should be 
noted that we only used this function to numerically analyze 
the color changes. Thus, the function might have no meaning 
in terms of the reaction mechanism. The obtained values of 
and RMSE for each model are summarized in Table 3. 
From the results obtained above, the Arrhenius plots 
were described using the autocatalytic model and the best- 
fit model. Figure 3 shows the autocatalytic model fitted to 
the changes in AE* ab at each temperature. High values 
of R ~ were obtained for all treatment temperatures 
(R- > 0.988). However, in spite of the high values of R , 
the experimental data deviated from the saturation of sig¬ 
moidal curve characteristic with the autocatalytic model as 
shown in Fig. 3, and the autocatalytic reaction is not likely 
to occur in the ventilated conditions as in the studies like 
ours [22, 28]. Using the rate constant k of each color 
parameter at each temperature calculated according to 
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Table 3 Comparison of coefficients of determination (R ) and root mean square errors (RMSE) of regression curves among several kinetic 
models and a best-fit model to changes in the color values during the treatment at 180°C 



Zero-order 

First-order 

Second-order 

Autocatalytic 

Best fit a 


R 2 

RMSE 

R 2 

RMSE 

R 2 

RMSE 

R 2 

RMSE 

R 2 

RMSE 

L* 

0.8846 

6.201 

0.9226 

5.174 

0.9533 

4.068 

0.9892 

1.901 

0.9999 

0.138 

< 3 * 

0.8665 

2.325 

0.7206 

3.190 

0.5660 

3.824 

0.9970 

0.351 

0.9998 

0.088 

b* 

0.5408 

7.816 

0.4369 

8.403 

0.3541 

8.837 

0.9924 

1.005 

0.9996 

0.226 

AE* ab 

0.7942 

9.987 

0.6683 

12.141 

0.5433 

13.782 

0.9878 

2.429 

0.9999 

0.205 


Underlined values are above 0.9900 in R 2 and below 1.000 in RMSE 
a Equation (14) was used as a best-fit model 



log (Treatment time) [h] 


Fig. 3 Change in AE* ab and its model: autocatalytic model fitted to 
each temperature 

Eq. (8), the activation energy E a was calculated from 
Eq. (9). The values of E a ranged from 101 to 125 kJ/mol, as 
shown in Table 4. Figure 4 shows the data superposed to 
the best-fit model using a T at each treatment temperature. 
The regression and the superposition fit the data set 
extremely well (R > 0.998). The activation energy E a was 
calculated according to Eq. (11) and the values of E a 
ranged from 118 to 125 kJ/mol. The values of E a for L* 
and AE* ab calculated from the autocatalytic model were 
approximately 20 kJ/mol smaller than those calculated 
from the best-fit model. This may be because the experi¬ 
mental data were a little outside the range of autocatalytic 
model under the present experimental conditions, as men¬ 
tioned above. Therefore, it is better to use the best-fit model 
with the TTSP method to accurately predict the color 
changes during heat treatment. 

Figure 5 shows the Arrhenius plots of AE* ab and their 
linear regression line as an example using the best-fit 
model. The regression lines for Arrhenius plots showed a 
good linearity for all color parameters (R > 0.987), sug¬ 
gesting that kinetic analysis can be applied to cellulose 
color properties and the prediction of color changes. The 
results also suggested that the color changes during heat 
treatment at the temperatures ranging from 90 to 180°C can 
be explained as the same reaction mechanism, presumably 


Table 4 Comparison of apparent activation energies (£ a ) calculated 
from the Arrhenius equation 



E a (kJ/mol) 


Autocatalytic 

Best fit a (TTSP) 

L* 

106 

125 

< 3 * 

125 

124 

b * 

113 

118 

AE* ab 

101 

120 


a Equation (14) was used as a best-fit model 



log (Treatment time) [h] 

Fig. 4 Change in AE* ab and its model: best-fit model and superposed 
data using a T at each temperature 

thermal oxidation, which occurs in the presence of oxygen 
and the absence of moisture. The values of E a calculated 
from the color parameters were similar to the reported 
values obtained from the other properties of cellulose, such 
as the degree of polymerization and mechanical properties 
[14, 21, 24, 29]. However, these reported values have wide 
variation and are still controversial [21]. 

Figure 6 shows the weight reduction as a function of 
time. The weight was reduced with increasing time at 120, 
150, and 180°C and slightly reduced at 90°C. As shown in 
the representative case in Fig. 7, the color changed 
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Fig. 5 Arrhenius plots of the change in AE* ab using the best-fit 
model (Eq. 14) and TTSP method 



log (Treatment time) [h] 

Fig. 6 Weight reduction during heat treatment as a function of time 

drastically before the remarkable weight reduction started, 
indicating that the color changes of cellulose are mainly 
due to the denaturation of cellulose and/or a slight amount 
of decomposition products of cellulose rather than the 
major decomposition of cellulose. The factors responsible 
for the color changes of cellulose were estimated to be 
carbonyl groups in the cellulose chains and the formation 
of colored low-molecular substances such as furan-type 
compounds [10, 30]. Lojewski et al. [9] reported that the 
increase of AE *ab had a significant linear relationship with 
the amount of furfural formed during heat treatment of 
cellulose at 90°C. Matsuoka et al. [31, 32] concluded that 
colored substances were formed at the reducing ends of 
cellulose during heat treatment at 60-280°C. These results 
suggested that the color changes were due to the denatur¬ 
ation of the reducing ends, the formation of furan-type 
compounds, or a combination of both these effects. The 
results obtained from this study will contribute to a better 
understanding of the degradation of cellulose and provide 



Weight reduction [%] 

Fig. 7 The relationship between AE* ab and weight reduction 

further evidence of the effectiveness of color parameters as 
an index of cellulose degradation. 

The contributions of cellulose to color changes of ligno- 
cellulosic materials were estimated by comparing them to 
color changes of hinoki (Chamaecyparis obtusa Endl.) 
wood which was treated and analyzed in almost the same 
conditions and procedures with the present study [33]. The 
values of L* decreased and < 2 *, /?* and AE* ab increased both 
in cellulose and wood. But in the later reaction stages of 
hinoki wood, the values of < 2 * and /?* started to decrease. 
The increase and subsequent decrease of < 2 * and /?* were 
also observed in the other wood species [34, 35]. The 
changes in L* and AE* ab of hinoki wood were similar 
degree to those of cellulose. These results show that cellu¬ 
lose is a major component that contributes to color changes 
of wood, however, other wood components such as 
extractives, lignin, hemicelluloses and minor components 
also have effects on color changes of wood, especially on 
the chromatic characteristics. The successful kinetic anal¬ 
ysis indicated that the thermal oxidation might cause the 
color changes of wood as well as cellulose, irrespective of 
treatment temperatures. The values of activation energy 
calculated from each color parameters of wood were lower 
than those of cellulose in all color parameters. The lower 
activation energy may be due to wood components other 
than cellulose. While color changes of ligno-cellulosic 
materials are due to the multiple and/or interactive effect of 
their components, the color changes of cellulose certainly 
contributes to that of ligno-cellulosic materials. 

Conclusions 

Color changes of cellulose heated at 90-180°C were 
measured and expressed by CIELAB color parameters. The 
value of L* decreased and the values of < 2 *, b* and AE* ab 
increased during heat treatment at all the temperatures. 
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Among the several kinetic models applied to the experi¬ 
mental data, the autocatalytic model best explained the 
changes in all color parameters. However, according to the 
fitting of the autocatalytic model and a review of the results 
in the literature, it is concluded that color changes during 
heating could not be explained by an autocatalytic reaction. 
To accurately predict the color changes during heat treat¬ 
ment, it was better to use the best-fit model with the TTSP 
method. The Arrhenius plot showed good linearity in the 
best-fit model, implying that the color changes might be 
induced by the same mechanism at the temperature range 
from 90 to 180°C. The apparent activation energies cal¬ 
culated from color parameters using the TTSP method 
ranged from 118 to 125 kJ/mol, which are similar to the 
reported values obtained from other chemical and 
mechanical properties of lingo-cellulose materials. 
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